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tions. Extensions for these models preserving this symmetry are studied using 



reproducing the known results in the weak coupling (band insulator) and 
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>-^j ■ strong coupling regime, first studied by Scalapino, Zhang and Hanke. Finally, 

a: 

the compact form of the variational wave functions allows to compute various 



correlation functions for these systems. 

I. INTRODUCTION 

The use of symmetries is an important tool to understand the effects of strong correla- 
tion in electronic systems. Recently, the S'0(3)-symmetry of the antiferromagnetic (AFM) 
order parameter has been combined with that of d-wave superconductivity to form a five- 
component vector order parameter!! It has been argued that the low energy sector of the 
resulting theory exhibits an approximate SO (5)-symmetry which allows to explain certain 
features such as the vicinity AFM order and superconductivity in the phase diagram of the 
high-T c materials. Numerical diagonalization studies have been performed and the spectrum 



of low lying excitations could in fact be classified according to this symmetry. 

A complementary approach has been the attempt to construct microscopic electronic 
systems with manifest SO (5) invariance and studies of such models to extract the low en- 
ergy behaviour. Scalapino, Zhang and Hanke succeeded in constructing a two-chain ladder 
Hamiltonian of this type and studied the strong coupling phase diagram of this system 
where they were able to identify several distinct phases (Ref. |^, referred to as SZH in the 
following). The properties of these systems at weak coupling in the metallic regime has 
been studied by means of bosonization&Q. Such ladder systems, in particular for magnetic 
insulators have attracted much attention recently due to the existence of various experi- 
mental realizations in materials closely related to the high-T c substances^. An interesting 
observation of Ref. is the existence of an SO (5) superspin phase which has been studied in 
a variational approach based on finitely correlated matrix product states similar to the ones 
used for S = 1 Haldane magnetsQia. Finitely correlated states have also been considered in 
electronic systems to describe aspects of the phase diagram of extended Hubbard modelsra!2l 
and other one-dimensional electronic models^JO. 

For SU(2) spin systems the variational approach has been generalized to lattices with 
ladder geometry and proven to give access to large parts of their phase diagramll^rta. This 
is the motivation for the present work where we extend the matrix product states originally 
introduced in Ref. ||to describe the strong coupling physics of the 50(5) superspin phase. We 
construct manifestly SO (5) invariant many particle wave functions from matrices containing 
all 16 electronic states on a given rung of the electronic ladder. The relative weight of 
the six different SO (5) multiplets on a rung is controlled by free parameters which are 
used to perform a variational study of the zero temperature phase diagram of the ladder 
at half filling. At strong coupling the results known from Ref. |2] are reproduced within our 
approach. Furthermore, at weak coupling and sufficiently large interchain hopping amplitude 
t± the matrix product state correctly describes the gapped ground state of a band insulator 
corresponding to a filled Fermi sea of electrons with one parity. For intermediate coupling 
we find a phase with finite amplitude of the 5*0(5) spinor quartets which are essential for 



the presence of a metallic phase of the ladder. The compact form of the variational states 
allows to study various correlation functions of interest. 

In the following section we present the classification of the electronic states of a two-leg 
ladder system according to the 50 (5)-symmetry and discuss all possible 5*0(5) symmetric 
single rung interactions. In Section [TIT] we review the SZH-model and consider tensor prod- 



ucts of rung states to include couplings of neighboring rungs. Section |V] deals with various 
SO (5) symmetric extensions of this model and a general construction routine for systems 
with exact finitely correlated ground states is given. Section [V| contains a detailed analysis of 
the ground state phase diagram of the system in the case of weak and intermediate coupling 
within a variational approach based on such wave functions. Furthermore we calculate the 
corresponding correlation functions within this approach. A summary of our results is given 
in Section IVTl 



II. ELECTRONIC STATES OF SO(5)-SYMMETRIC LADDER MODELS 

We consider a two-chain electronic ladder model with canonical creation and annihilation 
operators c\ (x), c a (x) for electrons (with spin-projection a =\, [) on sites x of the upper 
leg and analogous operators d\ (x) , d a (x) for the electrons on the lower leg . In order to 
discuss the 50(5) symmetry of the ladder model and to classify all the 16 possible states 
on a rung according to this symmetry, these operators are combined into four- dimensional 
50(5) spinorsi'0 

*a (z) = ( c T (x) , q (x) , 4 (x) , 4 (x) ) (a; even) (2.1) 

and 

** (s) = ( rf T (x) , d i (x) , c{ (x) , c\ (x) ) {x odd). (2.2) 

Using this definition the ten local generators L a b of the SO (5)-algebra on a single rung x 
are defined as 



L ab (x) = --%(x)T3p 



x 



a, b = 1, . . . , 5 



(2.3) 



Here T ab are ten antisymmetric, 4x4 matrices (their explicit form is given in Appendix |A|). 
A convenient basis of the Hilbert space on a single rung is diagonal in the quadratic Casimir 
charge 



C(x) = J2L 2 ab (x). 



(2.4) 



a<b 



In addition we choose to diagonalize the total charge Q = \{c^c + d)d — 2) and the z- 
component of the spin S z = \{$o z c + d)a z d). Based on the eigenvalues of C the Hilbert 
space can be decomposed into six 50(5) multiplets: 



Three SO (5) singlets (C = 0), for R see (TA^ 



hO = |n> ^ c MlAA | 0) 
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(2.5) 



An 5*0(5) vector quintet (C = 4) containing the ferromagnetically polarized state at 
half filling 
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(2.6) 



Two 50(5) spinor quartets (C = 5/2) for an odd number of electrons on a given rung 
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a = 1, ...4. 



(2.7) 



We label the states l^a) on a rung by the dimension d of the corresponding multiplet 
(a = 1, . . . , d) and an additional index k. Similarly, we can characterize product states on 
two rungs (see Sect, [i TTJ) . Alternatively, the vector quintet ( |2.6|) can be constructed from 
5-0(5) spinors |tfg (2) ) = ^(m ±n B )|fi), |*$ (4) ) = ^(n 2 ±n 3 )|fi) and |*$) = n 4 |0), 
with the superspin vector 



n a (x) = -*t ( X ) r^ ( x ) 



(2.8) 



Again, the explicit form of the 4x4 Dirac T-matrices T a is given in Appendix |A[ 

Any electronic ladder model with a local S'0(5)-symmetry on a rung has to preserve 
the degeneracy of the energy within the states of each single multiplet. The invariant 
Hamiltonian on a single rung can therefore be written as a sum over projection operators 
on these states: 

k = a 5 E i*£x*Si + E *i fc, ° E i*2x*2j + E A^l*?ix*& (2.9) 



M=l 



fc,J=l 



ti=l 



kl=l 



( (M) 



,('.*)> 



where A^' t; = (A^' )* because of the hermiticity of h x . All S'0(5)-symmetric terms on 
a rung can be expressed using linear combinations of these projection operators, e.g., the 
projection operator on the first singlet I^/qo) ls 



Hi = |tt$><*$l = -\s c (x)S d (x) + ^(S c (x)S d (x)) 2 , 



(2.10) 



with P d ' = fyd 1) {ipdLl anc ^ Sc(x) = \c^{x)ac{x). A complete classification of these terms 
is given in Appendix [B]. As a simple example we choose 
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which leads to the Hubbard-type Hamiltonian@ with an 5'0(5)-symmetry introduced by SZH 



-Hriing " Coulomb i -H Hopping 

= J2[ U (Kt(*) - ^){n ci {x) - -) + (c -+ d)) 

+ 7 K(x) - l)(n d (x) - 1) + J S c (x)S d (x) (2.11) 

-2t ± (4(x)rf CT (x) + .£f.c. 



where J = 4(U+V). This condition on the exchange amplitude guarantees the degeneracy 
between the states in the S"0 (5)-quintet and therefore the local S'0(5)-symnietry of the 
system. We will discuss this model and 50(5) symmetric extensions in the following sections. 

III. COUPLING OF NEIGHBORING RUNGS 

In order to describe an extended quasi-one dimensional electronic system one has to 
include coupling of neighboring rungs in addition to single rung interactions considered in 
the previous section. The simplest possible term is an 5*0(5) symmetric hopping term 
between adjacent rungs 

- 2*|| £ [4 (x)cv(y) + dt(x)d a (y) + H.c] , (3.1) 

(x,y) 

which can be brought into a manifestly S'0(5)-symmetric form using the alternating defini- 
tions of the spinors ( |2.1|J2.2j ). This hopping term together with the local rung interactions 
( P- 11| ) yield the complete SZH-modela. The ground state phase diagram of this system in 
the limit of strong coupling (U, V 3> t±, ty) has been determined by SZH using perturbation 
theory (see Fig. |2|). Four different phases have been established at half-filling: 
In phase I (occuring for < V < —2U) the model can be mapped onto an Ising-like system 
in a magnetic field: phase I a (V > —U/3) is a CDW phase and lb (V < —U/3) corresponds 
to the disordered Ising phase. Phase II is a spin-gap <i-wave phase (product of rung singlets), 
emerging for V > -U,U > and for V > -2U, U < 0. The phase III (V < -U, V < 0) 
is the superspin-phase where the 5*0(5)-quintet is dominant. For a further examination of 
this superspin phase, SZH have used the finitely correlated wave function 



^ s ZH ) = Ti(f[T a n a \n)] (3.2) 

(summation over the index a is implied and the trace is taken in the 4x4 matrix space 
where the T a are defined). In this form periodic boundary conditions have been imposed. By 
adding many particle interactions to their original Hubbard-type Hamiltonian this state fl3.2p 
can be made to be the exact ground state of the resulting model. This state has been argued 
to capture the essential physics of the superspin phase — similar to the role of the AKLT- 



model as a representative for a Haldane-gapped spin-1 chain. The wave function ( |3.2j ) will 
be the starting point for constructing a generalized matrix product wave function including 
all 16 states on a rung (see section |TV| ) and later be used for a variational study of the 
ground state phase diagram of the SZH model and its various SO (5) symmetric extensions 
beyond strong coupling (see section 0). The hopping term ( |3.1|) is one of many possibilities 
to include interactions between two adjacent rungs of the ladder but the requirement for 
a local 50 (5)-symmetry puts constraints on the explicit form of these terms. Explicit 
expressions for some of the interaction terms are listed in terms of electron operators in 
Appendix [B[ For a classification of these additional interactions we consider products of 
wave functions on two neighboring rungs x and y. A decomposition into 50 (5)-multiplets 
similar to (|2.5| - |2.7|) gives 50 different multiplets invariant under the action of the SO (5) 
generators £ a b(x, y) = L a b(x) + L a b(y). Tensor products containing a singlet factor on one of 
the rungs are trivial leading to simple product states, e.g. the 5*0(5) singlets |^oo)x|^oo)i/- 
Altogether there are nine singlets, 12 quartets and six quintets of this form. The remaining 
169 states are obtained by forming tensor products of quartets Q2.7]) and quintets (|2.6|). The 
decomposition of these products into irreducible representations of 5*0(5) reads 

4©4 = 1©5©10, 
4® 5 = 4©16, 

5©5 = 1©10©14 



(numbers denote the dimension of the corresponding SO (5) irrep). For example, one of four 
SO (5)-singlets in the tensor product of quartet states (|2.7|) is 



i*$W)> = \ (-i*2>i*8> - i*2>i*$> +i*S>i*S> + i*2>i*S>) • ( 3 - 3 ) 

Similar combinations of the rung states appear in the other states, the Casimir charges of 
the new multiplets are C = 6 for the decuplets, C — 10 for the 14-dimensional and C = 15/2 
for the 16-dimensional representations. The multiplets can be classified further according 
to the different eigenvalues of Q and S z on their member states. In Fig. [I] the state content 
of the various multiplets is shown. In the following we use this classification of the £0(5)- 
multiplets to construct ladder systems with exact ground states including different £0(5) 
symmetric nearest neighbor interactions. 

IV. EXTENSIONS OF SZH 

As mentioned in the Introduction the finitely correlated wave functions originally intro- 
duced to discuss the spin-liquid phases arising in one-dimensional higher spin Heisenberg 
models^ have recently been generalized to more general lattices. In particular, ladder mod- 
els whose ground states are of this form have been constructed^lLa'Ej. In these spin systems 
the ground state is of the form |\Po) = Ylx=i9x where g x is a (2 x 2) matrix containing the 
different states on a single site or rung x (e.g. spin-1 states for the AKLT model, singlet 
and triplet states for a two- leg S — 1/2 ladder, etc.). Different properties under translation 
in the extended direction can be realized by an appropriate choice of the free parameters 
appearing in g x (e.g. an alternation to introduce dimerizationll§E3El). Within a transfer 
matrix approach it is straightforward to compute various ground state correlation functions 
for different boundary conditions, periodic ones correspond to taking the trace of the matrix 
product wave function!^. 

For a further analysis of the SZH- model and the construction of £0(5)-symmetric ladder 
systems with exact ground states in matrix product form we have extended the wave function 
Q3.2|) to include the three SO( 5)-singlets (|2.5|) and the two SO (5) spinor quartets 

8 



i*o> =*& (n &({»})) ■ ( 4i ) 

Now (yf a . is a 5x5-matrix and p, (z = 1, ...,6) are variational parameters assigning different 
weights to the multiplets (|2.5| ) - (|2.7| ) on a rung (see Appendix [AJ). We restrict ourselves 
to the translational invariant case, where the parameters pi are chosen to be independent 
of the rung position x. In this case the matrix product wave function on two neighboring 
rungs contains two SO(5) singlets, two quartets, one quintet and one decuplet. The 14- 
dimensional and 16-dimensional representations are absent by construction. The states of 
the matrix product are linear combinations of the basis in Sect. [TiT] above, their explicit form 
is rather complicated. With respect to the spin-SU(2) subalgebra the remaining multiplets 
present in the matrix product contain spin singlet, doublet and triplet states only (states 
with total spin polarization S z > 1 are members of the 14 and 16-dimensional representation, 
see Fig. [1]). An immediate consequence is that the ansatz cannot be expected to describe 
the formation of ferromagnetic domains with higher spin states. An analogous argument 
holds for higher values of the charge Q, corresponding to strong local deviations from half 
filling. 

There is a simple way to construct spin ladder systems with matrix product wave func- 
tions as ground statescJ and a generalization to electronic ladder models with an SO (5)- 
symmetry is straightforward. The starting point is a general S'0(5)-symmetric Hamilton 
operator on two neighboring rungs 

4 16 14 

h x , x+ i = 2J -^16 E A^ + E ^ 14 Pu,n 

k,l=l £t=l £t=l 

5 10 10 5 

+ E Af '° E At + E 4 M) E P& (4-2) 

k,l=l £t=l k,l=l H=l 

16 4 14 

+ E a?'° E Pit + E * { o' l) ro% 

k,l=l /i=l k,l=l 

where P d ' = \wd,p) {^d,^] are projection operators on all possible S'0(5)-multiplets (see 
section [JT1]). The states \ip d l) are product wave functions on two rungs, k and / label 
the multiplet, // the states in the multiplet and d is the corresponding dimension of this 



irreducible representation. The hermiticity of h x , x +i requires A^ — (X d ' )* for the coupling 
constants, leaving altogether 322 free parameters in the Hamiltonian. A Hamiltonian H = 
J2 X hx,x+i has a finitely correlated ground state |^ ) — T\x9x with zero energy provided that 
the following conditions are satisfied^: 

• h XjX +x has to annihilate all states contained in the matrix elements of the product 
9x9x+\ 

• all other eigenstates of h XjX +i have positive energy 

Starting with an ansatz for g x in (|4.1| ) , one has to identify all multiplets \ipdj) contained 
in the product wave function g x g x+ i. These multiplets are labelled by indices id = 1, ■ ■■-.gd 
where the maximum number g^ is the number of multiplets with an equal Casimir-charge 
(d is the dimension of the irreducible representation), e.g. gio = 2 if the product wave 
function on two neighboring rungs contains two independent 5*0(5) decuplets. After the 
determination of the multiplet content of |\I/o) the corresponding parameters A^ ' in h XjX+1 
are set to zero to fulfill the first condition. The remaining operators in ( |4.2|) will now project 
on states not included in the matrix product wave function, which leads to zero energy for 
the ansatz. To satisfy now the second condition the reduced matrices A^ (I ^ id) have 
to be chosen positive definite (i.e. positive eigenvalues) such that ( [4.1| ) will be the lowest 
energy state of the system. 

In principle a general Hamiltonian where our ansatz (|4.1|) is the exact ground state, 
can be built by operators projecting on the other remaining S'0(5)-multiplets (e.g. the 
14-dimensional and the four 16- dimensional representations) and it has 249 free coupling 
constants A^ (/ 7^ id). We give explicit expressions for some of these operators in terms 
of electron operators in Appendix |B3j . In general however, the structure of these projection 
operators is quite complicated making it difficult to motivate these exactly solvable systems 
on physical grounds. 
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V. VARIATIONAL STUDIES OF THE PHASE DIAGRAM 

An examination of the SZH model beyond strong coupling can be done by using ( f4.1|) as 
a variational wave function. This wave function leads to the variational energy 

E rung = (V \H Coulomb \S> ) ~ (p 2 2 + pj)(U/2 -V) + 5pl(U/2 + V)- p\(Uj + 3V) (5.1) 

for the spin and charge interaction on a rung (see equation |2.12| ). Here p§ is the parameter 
corresponding to the 50(5)-quintet, p\ is the weight of the singlet \Q) and p 2 ^ of the 
symmetric and antisymmetric linear combinations of the other S'0(5)-singlets (|2.5| ). The 
variational energy corresponding to the hopping term on a rung is 



8 PlP2 + 2 {^f\ (w - h) 



(5.2) 



E t± = (^fo\HHopping\^o) ~ ^J 

and between two neighboring rungs (see equation |3~T|) it is 
E t {] ~ t\\ [p 4 p 5 (25pg + p\ - p\ - pfj + h 2 (p 2 Ps ~ 5pip 6 ) + (pl ~ pf) (5p 2 Pe - PiPs)} ■ (5.3) 



where hi = 5j9g + pf + p 2 + p\i h 2 = p\+ p\ and w = Jh\ + lQh 2 . 

In the strong coupling limit one can neglect these hopping terms ( |5.2j , [5.3| ). Minimizing 



the energy with respect to the Pi reproduces exactly the phase diagram calculated by SZH 
within perturbation theory (Fig. |^) where the phases are fixed by the largest amplitude Pi 
of the corresponding state and the crossover is continuous. The phase I is dominated by the 
bonding singlet state with amplitude p 2 , phase III is the superspin phase (p 6 ) and phase 
// consists of products of rung singlets (pi). In this approach with translationally invariant 
Pi the crossover between the two Ising-phases can not be reproduced. We now extend this 
analysis of the phase diagram of the SZH model to weak and intermediate coupling. 

A. Weak coupling phase diagram 

The band structure of the non-interacting system at half-filling is well known. For 
U= V= there are two energy bands, given by 

11 



e±(k) = ±2t ± - 4t|| cos(ife), -7r < k < n (5.4) 

and two different cases have to be distinguished: 

For t± < 2t|| the Fermi energy intersects the two bands (see Fig. ^.a) and for t± > 2tn 
they are separated by an energy gap (see Fig. |3|.b). The gapless system (t± < 2t\\) has been 
studied using bosonization of the low lying modes in the vicinity of the four Fermi points to 
obtain the phase diagram for weak coupling (U, V <C t±,tn): Lin et al.a found that at half 
filling the system is driven towards an integrable SO(8)-symmetric Gross-Neveu model in 
a weak coupling renormalization group analysis and predicted the occurrence of additional 
phases compared to the strong coupling case. The ongoing debate of these results (see the 
criticism of Ref. ^) cannot be clarified within the present ansatz: Using finitely correlated 
wave functions always leads to an exponential decay of correlation functions, indicating the 
existence of an energy gap between the ground state and the first excited state. 

B. Phase diagram for tj_ > 2tii 

For t± > 2t\\ the variational ansatz gives the exact ground state for the non-interacting 
system (U = V = 0): 

Choosing p\ = — -j^,P2 = 775 and Pi = for i = 3, ..., 6 we find 

L 

1^0) ~ n (r c \ d \ + c \ d \ - 44 - 44) (^ l°) ' ( 5 - 5 ) 

x=l 



which corresponds to complete filling of the modes with energy e_(fc) in (5.4), the band 



insulator. Consequently, we expect the variational approach to give reasonable results for 
the weak coupling phase diagram in this regime of hopping amplitudes. The quality of the 



approach can be measured by the mean deviation J {(AH) 2 ) = J (H 2 ) — (H) 2 of the energy. 
For U, V <C t±,t\\ the mean deviation stays small compared to the energy so that the ansatz 
should give reliable results. 

We find that only two phases are present in the weak coupling case (see Fig. |J) : the 
Ising phase / (j^) an d the spin-gap d-wave phase // (pi) already known from the strong 
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coupling diagram (see Fig. 0). The superspin phase disappears and also the 5'0(5)-quartets 
have no significant weight (pi,P5,Pe ~ 0) as expected for a band insulator. 

Considering the complete ground state phase diagram (see Fig. |j) we find an additional 
phase for intermediate coupling (U,V = t±,t\\) where the S'0(5)-quartets have the largest 
weight, in particular the rung-symmetric one \Q„) (W^)- Apart from these, the symmetric 
singlet state {^ a R a ^ ^ — \1/ a R a p^> p)\£l) -- which determines the ground state in phase I 
- has a significant weight. Due to the resonating structure of the ansatz and the relatively 
large variational value of {(AH) 2 ) the phase boundaries are not very accurate — for a more 
detailed study of this question the present work should be complemented by a numerical 
approach. As discussed earlier it is not possible within this approach to determine the 
position of the crossover line between the two Ising phases, or even whether this transition 
still occurs for the case of weak or intermediate coupling. 

C. Ground State Correlations 

The physics in the ground state is determined by ground state correlation functions, 
which are easily computed from matrix product wave functions. The matrix product ansatz 
( |4.1|) with the six free parameters Pi represents the ground state for a large class of models. 
We have calculated various correlation functions explicitly in the thermodynamic limit (L — > 
oo) for this variety of models (a detailed list is given in Appendix 0) and we determined the 
correlation length and the amplitude of different ground state correlations for the SZH-model 
when we used the ansatz as a variational wave function. 

The two-point correlations in matrix product states are always short-ranged (if not van- 
ishing) and have the following form 

(O\r)O(0)) = A({ Pi }) e"f 

They exhibit an exponential decay with the correlation length £ and amplitude A({p{\). As 
an example we consider the correlation length and amplitude of the expectation value of the 
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spin-spin correlation function (S Ct d(r)S Ct d(0)) (see Fig. ^) and of field correlators (c* a (r)cJ0)) , 
a,/3 G {#,«} (Fig. H) for the SZH-model on a circle in the [7-V-plane (with £7 2 + V 2 = 3) 
intersecting the phases /, II and the quartet phase (see Fig. f§). 

The spin-spin correlation function (S Cj d(r)S Ci d(Q)) in Fig. || is non-vanishing only in 
the quartet phase but with an extremely small correlation length indicating strong near- 
est neighbor correlations. For the electron-electron correlation in Fig. |6| with c 9iU (x) = 
(c|(x) + ci(x)) ± (<i|( x ) + ^l( x )) the correlation length £ is small for all angles but with a 
very large amplitude A except in the quartet phase. 

The sharp peak in both diagrams at ~ |-7r indicate the crossover of the phases II 
and I in Fig. ^ where the correlation length diverges. Calculating these correlations in the 
strong coupling limit the phase boundaries in Fig. ^ are denoted by very sharp peaks in 
the electron-electron correlation length £, + M m , with a non-vanishing amplitude. The 
spin-spin correlations are zero in the whole phase diagram and give no further hints of an 
underlying structure in the system. 

D. Variational examination of SO(5)-symmetric extensions 

Our variational approach is also suitable to study the phase diagrams of various SO(5)- 
symmetric extensions of the SZH model. We have considered additional interactions on a 
single rung and between two neighboring rungs, using the construction routine of section |II 
and|TV]. 

1. Single rung interactions 

All single rung interactions can be constructed using the projection operators of Sect. [I I 
and a detailed list of all possible terms can be found in the Appendix |||. Taking into account 
the operators P d ' with the coupling constants A^ ' leads to the following contributions to 
the variational energy (|5.1| ), calculated with the ansatz (|4.1|) 
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( P m) = f (5.6a) 

(P(ol } ) = ^(P2 + Ps) 2 (5.6b) 

(PiA)) = ^(P2 ~ Ps) 2 (5.6c) 



V2 



Wo) + Pfoin) = ^TPi(P3 + pa) (5.6d) 



».' 



(0,0) ^ - 1 (0,0) 






(0,0) -r ' (0,0)/ (( . 



E(Ai») = ^ i ^-P 5 ) 2 (8-88) 






2 
2 



BO = ^ • (mj) 

These modifications of the model cause some changes in the ground state phase diagram, 
e.g. the simple terms like (a) and (j) will only shift the phase boundaries without changing 
the general structure of the phase diagram. Other interactions like the pair hopping term 

(f) 

tpair (44 C T C I + h - C ) ~ ^(0,0) + ^(0,0) ( 5 - 7 ) 

will dramatically change the phase diagram (see Fig. [7]). 

For small negative values of t pa i r (\t pa i r \ < t\\) the Ising-phase I of the phase diagram 
in Fig. f| with the amplitude p 2 splits into two singlet phases: a symmetric (p 2 ) and an 
antisymmetric phase (^3) (see Fig. [7|.a , t pa i r = — 1), where the crossover line has the same 
gradient (U = — 2 V). Increasing the amplitude of t pair leads to a pure antisymmetric phase 
(see Fig. [7|.b , t pair = —4) in / and also to a strong change of the shape of the quartet phase. 
For small positive values of t pa i r the general structure of the phase diagram is preserved (like 
Fig. |J). In the regime of the coupling constants with t pair ^> t\\ the quartet phase vanishes 
(see Fig. ||). 
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Other interactions also exhibit strong effects on the phase diagram, e.g. including a 
quartet term like (i) which contains a hopping term on a rung and a bond-charge interaction 
(see App. FBI). 



tquar 



(c{d T + h.c.)(l-(n cl - n dl ) 2 )+ t<->J.l (5.8) 



leads to different phase diagrams, depending on the coupling constant. For positive values 
of t qU ar the quartet phase vanishes with increasing values of the coupling constant (like in 
Fig. D) until there are only the three known phases (see Fig. |^.a ). For t quar < the quartet 
phase grows (see Fig. ||.b ), dominated by the symmetric combination (p 4 ) of the states. 

The mean deviation in the weak coupling limit in these two special cases ( |5.7| , |5.8| ) is small 
compared to the energy (calculated on a circle with radius R = 0.1 around U — V — 0) 
except for the value t quar < — 1. The ansatz also provides very good results in the strong 
coupling limit (R > 100), except for the crossover lines to the superspin phase where the 
mean deviation is very large. The same problem occurs in the intermediate coupling regime 
in the quartet phase where the ansatz is not a good eigenstate of the system. 

We expect that including the other interactions on a rung will lead to similar changes in 
the ground state phase diagram. 

2. Two-rung interactions 

In most cases the SO (5)-symmetric interactions between two neighboring terms have a 
very complex structure but for some of them we can give simple expressions (see Appendix 
H). For them we can calculate the corresponding variational energy, e.g. the two-pair 
hopping term leading to an SO (5) singlet-singlet transition 



t 



2— pair 



dt(x)d\(y)d\(x)d\(y)cJx)cJy)c l (x)c l (y) + h.c (5.9) 



giving E two _ pair ~ 2t 2 - P air ($>i + Vz) (P3 - Vi) • The other S'0(5)-singlet interactions on two 
rungs lead to similar expressions, which will change the phase diagram (Fig. |]) in the Ising 
phase according to the value of the coupling constant t 2 - pa ir- 
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Including an SO (5)-quartet interaction in the SZH- model (see eqn. ( P18|) ) gives the 
variational energy 

^'" a^ + ftj fc-'*) (**-**)• (510) 

The phase diagrams obtained for different values of the coupling constant t qxy are very 
similar to the phase diagram in Fig. [|. The additional interaction has no significant effect 
except for minor changes of the crossover lines. 

VI. SUMMARY AND OUTLOOK 

We have constructed a large class of electronic ladder models with SO (5) symmetry 
having finitely correlated ground states and, consequently, correlation functions exhibiting 
exponential decay. These matrix product states have been used to perform a variational 
study of the ground state phase diagram of the SZH modelB for t± > 2t\\. For vanishing 
coupling the ground state of the band insulator is found to be in the class of variational states 
and at strong coupling the phases identified by SZH are reproduced. In the intermediate 
coupling regime signatures of a new phase dominated by local 50 (5)-quartets are found, and 
at weak coupling the 5*0(5) superspin phase is absent. Within our approach it is possible 
to compute various correlations giving further insights into the nature of the phases which 
have been identified. Finally, we have introduced various 5*0(5) symmetric extensions to 
the SZH model and discussed their impact on the phase diagram. In the future we will 
include dimerization in the matrix product ansatz for further studies of the Ising transitional 
in phase I and the possibility of spontaneous breaking of translational invariance in exactly 
solvable models. 
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APPENDIX A: THE GAMMA MATRICES AND THE VARIATIONAL WAVE 

FUNCTION 



For the construction of the SO (5)-invariant quantities, we have used the representation 
of the matrices in Ref. §. The five Dirac T-matrices have the following form 



/ 



r 1 



— ia, 



\ %o y 



\ 



/ 



r 



2,3.4 



(* \ 



a 1 



( 



r £ 



a, 



\ 



y<T y j 



(Al) 



where a are the Pauli matrices. The matrices T ab are defined by T ab = — |[r a , T b ] and the 
matrix R is given by 



R 



\ 



(A2) 



/ 



11 

-11 

Using these definitions, it is simple to construct the matrix g x of the variational wave function 
(|4. 1|) . It has the following structure 



/ 



/ 



\ 



(*). 



p 6 r a n a \n) + EUPi\n!o)*- 



\ 



la) 

I Q2) 

1 93) 

1 94) 



(A3) 



/ 



\-\qs) - I94) \qi) |<?2> 
The three SO(5)-singlets are included in this ansatz only on the main diagonal elements. 
^0,0) — l^oi) from Q2.5|) and l^oo ) are the symmetric and an antisymmetric combinations 
of the two other singlets 



Ko 3) > ~ (*^*J T VaRafiVfl) |0>. 



(A4) 



The quartets enter the matrix g x in \q a ) = P4\Qt)+P5\Q a ) where \Qt) are the symmetric 
and antisymmetric combinations of ( |2.7| ) 



Qa 



I 



± 



I 



T 



± 



T 



r 

Tl, 



± 



U' 
t, 



r 
u. 



± 



Tl' 
I. 



(A5) 



They are arranged in the right column and the lowest row of ( |A3|) in such a way that in the 
product g x g x+ i one has SO(5)-singlets on the diagonal only. 



APPENDIX B: SO(5)-SYMMETRIC OPERATORS ON ONE RESP. TWO RUNGS 

We present here a selection of various SO (5)-symmetric terms on a single and on two 
rungs ( [B 1| , |B 2| ) . Furthermore, a list of terms is given for which our matrix product ansatz 
Q4.1| ) would be the lowest energy state 



1. Single rung interactions 

We now present all possible S'0(5)-symmetric terms on a rung. Their general construc- 
tion is done in terms of projection operators on the different S'0(5)-multiplets. Expressed 
through electronic operators most of them are already known from the SZH-model (|2.11|) 
and an additional biquadratic exchange. As a shorthand notation we introduce 

[U,V,J,a]= U ^{x)- l -){n cl {x)- l -) + {c^d)^ 

+ V (n c (x) - l)(n d (x) - 1) + J S c (x)S d (x) (Bl) 

+a(S c (x)S d (x)) 2 

In addition we find various single electron and pair hopping terms together with bond-charge 
type interactions. Using the notation of Sect. || for the projection operators on a rung we 
obtain the following terms by projection on the singlets 

(B2) 
(B3) 
(B4) 
(B5) 

(B6) 
(B7) 



pl,l _ 

^0,0 — 

p2,2 _ 
-'0,0 _ 

p3,3 _ 
-'0,0 _ 

-^0,0 + J 

-Po,b + J 

pS:S + - 


^o,o/\^o,ol - 
■1 12 4" 

2' - 4'3'3. 
■1 12 4" 


r 1 4 " 


l-n c ) 

l-n c ) 

nci{ridi 


— c <-> d] 


L2' 4' 3' 3j 

p2,l _ 1 \( 

°>° " V2 t { 

p3,l _ ' 

p o,b = "t i C T 


2 L n rfT n dH 
ddr + h.c.)n c 

(cldt + h.c.) 
Cj_ + h.c. . 


1)+ W 
-i)+T<-» 


l" 



19 



The projection operators on the quartet states read 



ji,i 



0,0, 
0,0, 



8 16 



3' 3 

8 16 

3'~y 



+ (1 - n c] n ci )n d + n d] n di (n c - 2) 
+ (1 - n d ^n di )n c + n c ^n ci (n d - 2) 



4 
ZvM,/* 

E 3S + p ii = [( c R+ ^c)(i-Ki- ^) 2 )+ rw 



and finally, projection on the quintet gives 



E-^ 



1 13 20 

'2'T'y 



(B8) 

(B9) 

(B10) 



(Bll) 



2. Interactions between neighboring rungs 

Equivalently, the 50 (5)-symmetric expressions on two rungs can be classified. The 
choice of the basis on the two- rung system is very important for the structure of the 50(5)- 
symmetric terms. Using the simplest combination the product of an SO (5)-singlet on one 
rung and another 50(5) multiplet on the other gives for a projection operator e.g. 






(B12) 



k,l 



for the product of an 50(5)-singlet on rung x and an 50(5)-quintet on y, where P d is 



defined in section |T[ The numbers k and I in P d ' depend on the way the different multiplets 
on the rungs are labelled. Another example is an operator projecting on an SO (5)-singlet 
on each of the rungs 



,(2) 



(2), 



>1,1; 



,2.2, 



P$S(x,y) = \*$(x,y)){% d i(x,y)\ = P${z)P$(y). 



(B13) 

(see equation ( |2.5| ) for the definition of the wave functions). All projection operators of 
states consisting of at least one SO (5)-singlet on a rung can be decomposed in the same 
manner. For some of these operators a compact representation in terms of electron operators 
is possible. As an example consider the operator 
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Jf~\J 



x)c\(x)d\(y)d\(y)c^(y)c l (y)d^(x)d l (x) + h.c. 



(B14) 



describing pair exchange between two neighboring rungs. It causes a transition between two 
5*0(5) singlet states and can be written as 



P$X(v)P§S(x) + h.c 



Other SO (5) singlet-singlet transitions of this type are 



(B15) 



d\(x)d\(y)d\(x)d\(y)c^x)c^(y)c l (x)c i (y) + h.c 
N d(y)n c ^(y)n cl (y) [d\(x)d\(x)c^(x)c i (x) + h.c 
N c (x)n d ^{x)n di {x) \c\{y)c\{y)d^{y)d l {y) + h.c 



(B16) 



where 



N Q (y) = (1 - n a1 (y) - n al (y) + n a ^(y)n ai (y) 



a G {c, d} 



(B17) 



Similar terms are obtained from projection operators on direct products of an 50(5)- 
singlet on one and an 50 (5)-quartet on the other rung, e.g. 



(KtO) - n d](x)f - l) c\(y)c\(x)c\(y)d 1 (y)d i (x)d l (y) 
+ (WO) - n di(x)) 2 - l) c\(x)c\(y)c\(y)d 1 (x)d 1 (y)d i (y)] + h.c 



(B18) 



The projection operators on the remaining 169 states with a structure similar to ( |3.3|) 
cannot easily be decomposed in this way. They are significantly more complex, generi- 
cally their expansion into electronic operators produces complicated bond-charge interaction 
terms. Still, forming suitable linear combinations of such terms can lead to simpler 50(5)- 
symmetric terms on two rungs, e.g. the pair hopping term in ( p.l|) or a diagonal hopping 
term 



H [ d U x ) c °(y) - c l(x)d a {y) + h.c. 



(B19) 
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3. SO(5)-symmetric Hamiltonians with exact Ground States 

At the end of section [IV] we claimed that a general Hamiltonian where our ansatz ( |1.1| ) 
is the ground state of the system is given by 

4 16 14 

h x , x +i = Yl ^16 Y Pilv + Y ^14 Pu,n + additional terms. (B20) 

k,l=l /i=l /i=l 

The coupling constants have to be chosen such that A i4 > and the matrix Ai 6 of coupling 
constants is positive definite. This implies that E = is a lower bound on the spectrum and 



therefore the state (4.1) -- having zero energy by construction -- will be a ground state. 
The additional terms in ( |B20| ) are the projection operators on the remaining multiplets not 
present in the matrix product wave function. For example the projection operator on one 
of the 5'0(5)-singlets not present in this product reads 

-c\(x)c\(y)c\(x)c\(y)d^x)d^(y)d l (x)d l (y) + n d1 (x)n d1 (y)n dl (x)n dl (y)N c (x)N c (y) (B21) 
+ (c «-> d)\ 






Just as Ai6 above the matrices A^ (d = 0,4,5, 10), coupling the projection operators on the 
remaining multiplets have to be chosen to be positive definite. 

APPENDIX C: CORRELATION FUNCTIONS 

The calculation of expectation values between matrix product states is straightforward 



using a transfer matrix method (see e.g. Ref. p2|) : 



To this end we define a 25 x 25 transfer matrix G on a rung 

G ai: a 2 ~ G (h,ji),(i2,J2) — 9( iui2 )9(j u j 2 ) ( C1 ) 

with the indices 

«! = !,... ,25^ (11),. ..,(15), (21),. ..,(55). 

In terms of G the norm of the ground state can be written as 
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25 



<v]>o | vDo) = Tr G L = $>^ 



(C2) 



*=1 



where Aj are the eigenvalues of G. In the thermodynamic limit (L — > oo) the largest 
eigenvalue Ai dominates this expression and we obtain (vI/qI^o) ~ Af. Similarly, one-point 
correlators of an operator O are 



(0) = -( ei |Z(0)| ei ) 



(C3) 



and a two-point correlation function reads 



25 1 / X \ r 



)|e n )(e n |Z(O r )|ei). 



(C4) 



Here \e n ) are the eigenvectors with eigenvalue A n of G and Z(Oj) ~ g^Oig is the transfer 
matrix related to the operator Oj. With the matrix ( |A3| ) the largest eigenvalue of G is given 
by 



Ai = - {h + w) 



(C5) 



with /i x = 5pg + Pi + pi + Pi; ^2 = v\ + P5 an d w — \Jh\ + 16/z|. 

This enables us to calculate the expectation values of any operator acting on a single or 
two rungs, respectively. For example, we find 

w - hi + 8j9g 



(S) = 0, 
for local magnetic moments and 

(4(x)c /3 (x)) 



(S 



i\2 



Aw 



(C6) 



w — h\ , , 9 . 

^2 
w — hi 



h? 



(2p 4 Ps) - 4p 2 P3 



a^/3 



(C7) 

(C8) 



for electronic expectation values. Here, a,/?G {g, u] and c 9>u (a;) = (c^(x) + q(x)) ± (rft( x ) + 
di(x)). 

Correlations between the total spin on two rungs decay exponentially 



(SlS r ) 



'hi-4plY fw - hi + 8pf 



Aw(hi + w) \ Ai 



hi - 4p§ 



(C9) 
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as expected for finitely correlated states. Spin-spin correlations between individual sites on 
rungs separated by a distance r can be expressed as 

(S a (r)S,(0)} = A a ,({ Pl }) (^=^) + B a ,({ Pl }) (^=^) , (CIO) 

where the amplitudes A a/3 ({pi}) and B a/3 ({pi}) depend on the choice of a and (3, i.e. whether 
correlators of spins on the same or on different legs of the ladder are considered. Analogously, 
one can study electronic correlations, e.g. 

(cl(r)cl(r)c g (0)c u (0)) = -l&Sr) 

(cl u (r)c g:U (0)) = C g , u {{ Pl }) (j^J + D g , u ({ Pi }) (-^j . (Cll) 

with the amplitudes C 9)U and D gu . 
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FIG. 1. The irreducible £0(5) representations appearing on a pair of rungs decomposed cor- 
responding to the eigenvalues of Q and S z : (a) the quartet (with Casimir charge C = 5/2), (b) 
the quintet (C = 4), (c) the ten-dimensional (C = 6), (d) the 14-dimensional (C = 10) and (e) 
16-dimensional {C = 15/2) irrep (double circle indicate two states with identical eigenvalues). 
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V=-2U 
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III 



v=-u 



FIG. 2. Strong coupling phase diagram, U and V measured in units of t\\ 





(b) 
FIG. 3. Band structure of a two-leg ladder model for (a) t± = t\\ and (b) t± = 2t\\ 
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V=-2U 



V=-U/3 




V=-U 

FIG. 4. Phase diagram for t± = 2t\\ (t\\ = 1): the phase boundaries were calculated by com- 
paring the amplitudes of the different multiplets. 




FIG. 5. Correlation length and amplitude for the spin-spin correlation function , the full line 
corresponds to (S c(d) (r)S c ( d )(0)) and the dotted to (S c ^(r)S d ^)(0)). 
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FIG. 6. The correlation length and amplitude of the expectation value (ct u (r)c (0)), the full 
line corresponds to (cj, (u) (r)c g(u) (0)) and the dotted to (c* (u) (r)c u(ff) (0)). 
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FIG. 7. Phase diagram for t± 



v=-u (b) v=-u 

2tj| including pair hopping for (a) t pa i r = —t\\ and (b) 



tpair = — 4tii with in = 1. 
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FIG. 8. Phase diagram for t± = 2t\\ including pair hopping for t va ir = 4tn with t\\ = 1. 
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FIG. 9. Phase diagram for t± = 2t\\ including a quartet term with (a) t qucir = +4iy and (b) 




tquar = — 4in with t» = 1. 
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